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Abstract 

For years limited Monte Carlo simulations have led to the suspicion that the 
time to failure of hierarchically organized load-transfer models of fracture is 
non-zero for sets of infinite size. This fact could have a profound significance 
in engineering practice and also in geophysics. Here, we develop an exact 
algebraic iterative method to compute the successive time intervals for indi- 
vidual breaking in systems of height n in terms of the information calculated 
in the previous height n — 1. As a byproduct of this method, rigorous lower 
and higher bounds for the time to failure of very large systems are easily 
obtained. The asymptotic behavior of the resulting lower bound leads to the 

evidence that the above mentioned suspicion is actually true. 
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There are few mechanical problems more complex and difficult to cast into a definite 
physical and theoretical treatment than the range of phenomena associated with fracture. 
Furthermore, there are few problems with a wider field of application: material science, 
engineering, rock mechanics, seismology and earthquake occurrence. Our understanding of 
fracture processes in heterogeneous materials has recently improved with the development of 
simple deterministic and stochastic algorithms to simulate the process of quasistatic loading 
and static fatigue Q. The load-transfer models used here are called fiber-bundle models, 
because they arose in close connection with the strength of bundles of textile fibers. Since 
Daniels' and Coleman's |§ seminal works, there has been a long tradition in the use of these 
models to analyze failure of heterogeneous materials ||. Fiber-bundle models differ by how 
the load carried by a failed element is distributed among the surviving elements in the set. In 
the simplest scheme, called ELS for equal load sharing, the load supported by failed elements 
is shared equally among all surviving elements. Important from the geophysical point of view 
is the hierarchical load-sharing (HLS) rule introduced by Turcotte and collaborators in the 
seismological literature ||. In this load-transfer modality the scale invariance of the fracture 
process is directly taken into account by means of a scheme of load transfer following the 
branches of a fractal (Cayley) tree with a fixed coordination number c. The HLS geometry 
nicely simulates the Green's function associated with the elastic redistribution of stress 
adjacent to a rupture. In the static case, the strength of an HLS system tends to zero 
as the size of the system, N, approaches infinity, though very slowly, as l/log{logN) 0. 
The dynamic HLS model was introduced in the geophysical literature in reference ||; their 
specific aim was to find out if the chain of partial failure events preceding the total failure of 
the set resemble a log-periodic sequence 0. In the analysis of 0, it appeared that, contrary 
to the static model, the dynamic HLS model seemed to have a non-zero time-to-failure as the 
size of the system tends to infinity. In this paper we provide evidence that this conjecture 
is true by means of an exact algebraic iterative method where trees of n levels, N = c n , are 
solved using the information acquired in the previous calculation of trees of n — 1 levels. 
As a byproduct of this method, we obtain a rigorous lower bound for the time to failure of 
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infinite size sets which results in being non-zero. 

In H we showed how in fiber-bundle dynamical models, using a Weibull exponential 
shape function and a power-law breaking rule, one can devise a probabilistic approach which 
is equivalent to the usual approach || where the random lifetimes are fixed at the beginning 
and the process evolves deterministically. Without loss of generality, from the probabilistic 
perspective the set of elements with initial individual loads a = uq = 1 and suffering 
successive casualties is equivalent to an inhomogeneous sample of radioactive nuclei each 
with a decay width a p ; p being the so called Weibull index, which in materials science adopts 
typical values between 2 and 10. As time passes and failures occur, loads are transferred, a 
changes and the decay width of the surviving elements grow. In the probabilistic approach, 
in each time step defined as: 
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one element of the sample decays. The index j runs along all the surviving elements. The 
probability of one specific element, m, to fail is p m = u^S. Eq. ([I]) is the ordinary link 
between the mean time interval for one element to decay in a radioactive sample and the 
total decay width of the sample as defined above. Due to this analogy, radioactivity terms 
will frequently be used in this communication. Note that we use loads and times without 
dimensions. The time to failure, T, of a set is the sum of the iV 5s. For the ELS case, the 
value of 5 as defined in Eq. ([]]) easily leads to T = 1/p, which is the correct result for the 
time to failure of an infinite ELS set 0. In this communication we will apply these ideas to 
the HLS case obtaining the 5s algebraically without having to use Monte Carlo simulations. 

To give a perspective of what is going on in the rupture process of a hierarchical set we 
have drawn in Fig. [I] the three smallest cases for trees of coordination c = 2. Denoting by n 
the number of levels, or height of the tree, i.e. iV = 2™, we have considered n — 0, 1 and 2. 
The integers within parenthesis (r) account for the number of failures existing in the tree. 
When there are several non-equivalent configurations corresponding to a given r, they are 
labeled as (r, s), i.e, we add a new index s. The total load is conserved except at the end, 
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when the tree collapses. Referring to the high symmetry of loaded fractal trees, note that 
each of the configurations explicitly drawn in Fig. [I] represents all those that can be brought 
to coincidence by the permutation of two legs joined at an apex, at any level in the height 
hierarchy. Hence we call them non-equivalent configurations or merely configurations. In 
general, each configuration (r, s) is characterized by its probability p(r, s), YjP{ t i s ) — 1> an d 

s 

its decay width T(r, s). The time step for one-element breaking at the stage r, is given by 

This is the necessary generalization of Eq. (|l|) due to the appearance of non-equivalent 
configurations for the same r during the decay process of the tree. In cases of branching, 
the probability that a configuration chooses a specific direction is equal to the ratio between 
the partial decay width in that direction and the total width of the parent configuration. 
We will compute at a glance the 5s of Fig. [I] in order to later analyze the general case. To 
be specific, we will always use p — 2. For n = 0, we have T(0) = l 2 and do — -k — 1 = T. 
For n = 1, T(0) = l 2 + l 2 = 2, 5 = ±, T(l) = 2 2 , 5 X = \ and hence T = § + \ = §. For 
n = 2, T(0) = l 2 + l 2 + l 2 + l 2 = 4, 5 = J, T(l) = 2 2 + l 2 + l 2 = 6, 5 1 = |; now we 
face a branching, the probability of the transition (1) — > (2, 1) is I and the probability of 
the transition (1) -» (2,2) is |, on the other hand T(2, 1) = 2 2 + 2 2 = 8 = T(2,2), hence 
<^2 = |'| + |'| = |- Finally 5 3 = ^ and the addition of 5s gives T — ||. 

Now we define the replica of a configuration belonging to a given n, as the same config- 
uration but with the loads doubled (this is because we are using c = 2). The replica of a 
given configuration will be recognized by a prime sign. In other words (r, s)' is the replica 
of (r,s). Any decay width, partial or total, related to (r,s)' is automatically obtained by 
multiplying the corresponding value of (r, s) by the common factor c p = 2 P = 4. This is a 
consequence of the power-law breaking rule assumed. The need to define the replicas stems 
from the observation that any configuration appearing in a stage of breaking r of a given n, 
can be built as the juxtaposition of two configurations of the level n — 1, including also the 
replicas of the level n — 1 as ingredients of the game. In Fig. [I], one can observe the explicit 
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structure of the configurations of N = 4 as a juxtaposition of those of iV = 2 and its repli- 
cas. From this perspective we notice that as the configurations for the height n are formed 
by juxtaposing two configurations of the already solved height n — 1 and its replicas, with 
the restriction that the total load must be equal to 2 n because the total load is conserved, 
the single-element breaking transitions occurring can only be of three types. One case a) 
corresponds to the breaking of one element in a half of the tree while the other half remains 
as an unaffected spectator. Another case b) corresponds to the decay of the last surviving 
element in a half of the tree, which provokes its collapse and the corresponding doubling of 
the load borne by the other half. Finally, the third case c) corresponds to the scenario in 
which one half of the tree has already collapsed and in the other half one break occurs. In 
this third case the decay width is obtained from the information of the replicas. This holds 
for any height n, allows the computation of all the partial decay widths in a tree of height 
n from those obtained in the height n — 1 and will be illustrated in Fig. |3|. 

Using henceforth a convenient symbolic notation, in Fig. 0a we have built what we call the 
primary width diagram for n = 3 resulting from the juxtapositions of pairs of configurations 
of the previous diagram of n = 2 and its replicas. It is implicitly understood that time 
flows downward along with the breaking process. In Fig. |2|a the boxes represent n = 3 
configurations, and at their right, the two quantities in parentheses indicate the two n = 2 
configurations forming them. The value of the partial width connecting a parent and a 
daughter is written at the end of the corresponding arrow. The sum of all the partial widths 
of a parent configuration in a branching is always equal to the total decay width, T, of 
the parent. From this primary width diagram one deduces the probability of any primary 
configuration at any stage r of breaking, and consequently 5 r is obtained using Eq. 
Finally, by adding all the 5s we calculate T(n = 3). To illustrate the connection between 
the explicit configurations as those of Fig. [1| and the somewhat hermetic notation of Fig. |2|, 
in Fig. 0, we have shown three explicit examples relating them. Fig. |3] is selfexplanatory. 
The three examples correspond to the ), b) and c) mentioned before. 

By iterating this procedure, that is by forming the primary diagram of the n+1 height by 
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juxtaposing the configurations of the primary diagram of the height n, we can, in principle, 
exactly obtain the total time to failure of trees of successively doubled size. Two examples 

are T(n = 3) = and T(n = 4) = SotSS • The P roblem arises from the vast 

amount of configurations one has to handle in the successive primary diagrams. This fact 
eventually blocks the possibility of obtaining exact results for trees high enough as to be able 
to gauge the asymptotic behavior of T in HLS sets. That is why, taking advantage of the 
general perspective gained with the exact algebraic method, henceforth our more modest 
goal will be to set rigorous bounds for T. This task is much simpler. 

To set bounds, we will define effective diagrams in which for each r there is only one con- 
figuration which results from fusing in some specific appropriate way all the configurations 
labeled by the different s values; see Fig. ||b. These effective configurations are connected 
by effective decay widths denoted by a r . Thus, the effective diagram for any n is calculated 
from its primary diagram, and then is used to build a primary diagram of the next height 
n + l. For n — 0, 1 and 2, a r = T(r) because Vr, the T(r, s) are equal to a unique value T(r). 
The point is how to define a r for n > 3, so that all the 8(r)(n > 4) and hence T(n > 4) 
are systematically lower (or higher) than its exact result. This goal is easily accomplished 
using, 

a r = T max (r) (or T min (r)) (3) 

i.e., by assuming that the configuration of maximum (minimum) width saturates the single 
element decay of the stage r. Better rigorous boundings are obtained using the geometric 
mean (or the harmonic mean), namely 

« r =n r (^r M ^ = , \ i ) (4) 

The bounds obtained from these formulae are plotted in Fig. It is clear that those based 
on Eq. (|3p are poor, in fact the lower bound goes quickly to zero. On the other hand, those 
based on Eq. (f|) are stringent and useful. The value obtained for the lower (higher) bound 
to T, from Eq. (Q), will be called TJ (T h ). The arithmetic mean, i.e., a r = J2p( r , s ) ' T( r ; s ) 

s 
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also leads to a lower bound but it is not as good as that coming from the geometric mean. 
The bounds emerging from these three means are rigorous because when configurations of 
different decay width are fused, the form of the generated function appearing in the calculus 
of the 5s are concave (convex). This will be explained elsewhere. 

We have fitted the data points of T t by an exponential function of the form 7} = TJ j00 + 
ae~ b ( n ~ no \ T} j00 , a, b and no being fitting parameters. The success of this fitting is crucial 
because this exponential decay to a non-zero limit is the hallmark of our claim. By choosing 
subsets formed only by points representing big systems we observe a clear saturation of 
the asymptotic time-to-failure towards T/ j00 = 0.32537 ± 0.00001. An analogous analysis 
of T h leads to T Koo = 0.33984 ± 0.00001. Similar fittings of the Monte Carlo data points 
are inconclusive, due to the intrinsic noisiness of the MC results and the limited size of 
the simulated sets (N < 2 16 elements). What this result implies is that a system with a 
hierarchical scheme of load transfer and a power-law breaking rule (c = 2, p = 2) has a 
time-to-failure for sets of infinite size, T^, such that 0.32537 < < 0.33984. Thus, there 
is an associated zero-probability of failing for T < and a probability equal to one of 
failing for T > T^. The critical point behavior is thus confirmed. Invoking conventional 
universality arguments one deduces that this property holds for a hierarchical structure of 
any coordination. The case of dynamical HLS sets using an exponential breaking rule will 
be reported shortly. 
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P. thanks Diego Vazquez-Prada for discussions. Y.M thanks the AECI for financial support. 
This work was supported in part by the Spanish DGICYT. 



7 



REFERENCES 



* On leave from Departamento de Ffsica, Technological University of Havana, ISPJAE, 
Havana 19390, Cuba. 

[1] P.L. Leath and P.M. Duxbury. Phys. Rev. B49, 14905 (1994); H.J. Herrman and S. 
Roux (eds). Statistical Models for the Fracture of Disordered Media (North Holland, 
Amsterdam, 1990), and references therein. 

[2] H.E. Daniels, Proc. R. Soc. London A183, 404 (1945); B.D. Coleman, J. Appl. Phys. 
28, 1058 (1957). 

[3] D.G. Harlow and S.L. Phoenix, J. Composite Mater. 12, 195 (1978); R.L. Smith and 
S.L. Phoenix, J. Appl. Mech. 48, 75 (1981); D.G. Harlow and S.L. Phoenix, J. Mech. 
Phys. Solids. 39, 173 (1991). 

[4] R.F. Smalley, D.L. Turcotte and S.A. Sofia, J. Geophys.Res. 90, 1894 (1985). 

[5] W.l. Newman and A.M. Gabrielov, int. J. Fracture 50, 1 (1991). 

[6] W.L Newman, A.M. Gabrielov, T.A. Durand, S.L. Phoenix and D.L. Turcotte, Physica 
D 77, 200 (1994); W.L Newman, D.L. Turcotte and A.M. Gabrielov, Phys. Rev. E 52, 
4827 (1995). 

[7] D.Sornette and C.G. Sammis, J. Phys. France 5, 607 (1995); H. Saleur, C.G. Sammis 
and D. Sornette, J. Geophys. Res. 101, 17661 (1996). 

[8] J.B. Gomez, Y. Moreno and A.F. Pacheco, Phys. Rev. E, 58, 1528 (1998). 



8 



FIGURES 

FIG. 1. Breaking process for the three smallest trees of coordination c = 2 (N = 1, 2, 4). The 
integers in parenthesis (r) represent the number of breakings occurred. The 5s stand for the time 
steps elapsed between successive individual breakings and the numbers under the legs indicate the 
load they bear. 

FIG. 2. Width diagrams for the breaking process of N = 8. a) Primary diagram. The various 
configurations are labeled by the integers inside the boxes. The integers in parentheses at the right 
of the boxes represent the two N = 4 configurations juxtaposed to form each of the N = 8 config- 
urations. The number accompanying an arrow connecting two boxes stands for the dimensionless 
partial decay width of that transition, b) Effective diagram. This is obtained from a), and is used 
to obtain rigorous bounds for the next height, N = 16. The as are explained in the text. 

FIG. 3. Relation between the explicit and the symbolic notations for configurations of N = 8, 
built juxtaposing pairs of N = 4. See the text for details. 

FIG. 4. Dimensionless lifetime, T, for a fractal tree of height n: the circles are obtained from 
Monte Carlo simulations; the continuous lines are the stringent bounds, and the dotted lines are 
the gross bounds. 
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